T-Duality and the Spectrum of Gravitational Waves 
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In the inflationary universe scenario, the physical wavelength of cosmological fluctuation modes 
which are currently probed in observations was shorter than the Hubble radius, and in fact shorter 
than the Planck and string lengths, at the beginning of the period of inflation. Thus, during the 
early stages of evolution, the fluctuations are subject to Planck scale physics. In the context of an 
inflationary cosmological background, we examine the signatures of a specific modified dispersion 
relation motivated by the T-duality symmetry of string theory on the power spectrum of gravita- 
tional waves. The modified dispersion relation is extracted from the asymptotic limit of the string 
center of mass propagator. 

PACS numbers: 



I. INTRODUCTION 

The inflationary universe acts as a microscope which 
allows us, by means of cosmological observations per- 
formed today, to probe scales which at the beginning 
of the period of inflation were sub-Hubble (i.e. the 
wavelength was smaller than the Hubble radius), and in 
fact sub-Planckian (wavelength smaller than the Planck 
scale) 0. This is possible since during the period of 
inflation (which we for concreteness take to be almost 
exponential) the physical wavelength of fixed perturba- 
tion modes increases exponentially, whereas the Hubble 
radius remains approximately constant. Thus, in prin- 
ciple, it is possible to probe Planck-scale or string-scale 
physics using current observations. 

Crucial to the study of trans-Planckian effects on cur- 
rent cosmological observations is the ability to track 
the cosmological fluctuations from the earliest times 
and smallest scales to late times. Signatures of trans- 
Planckian physics then emerge via encoding the new 
physics which is effective at the Planck scale on the initial 
conditions and early evolution of the linearized fluctua- 
tions. This program was initially carried out 2, 3] mak- 
ing use of ad-hoc modified dispersion relations for the 
cosmological perturbations. These dispersion relations 
are linear (i.e. the usual ones) for physical wavenumbers 
smaller than a cutoff value which corresponds to the en- 
ergy scale of the new physics, and deviates for larger 
wavenumbers. This approach has the advantage that it 
is easily computationally tractable. It was shown '2*1 that 
observable effects are possible provided that the disper- 
sion relation violates the adiabaticity condition, i.e. that 
the early evolution of the fluctuations is non-adiabatic 
from the point of view of the evolution which would have 
been obtained had the dispersion relation been the un- 
modifled one. There are some constraints on the mag- 



nitude of trans-Planckian effects which come from back- 
reaction analyses |^], but these are less severe than ini- 
tially conjectured in j^Q. 

There are other approaches to studying the possible 
imprints of trans-Planckian physics on cosmological ob- 
servations. In a "minimal" approach, initial conditions 
on the fluctuations are imposed not at some early initial 
time which is the same for all wavenumbers, but for each 
wavenumber at the time when the corresponding wave- 
length equals the scale of the new physics 0. There are 
approaches which attempt to include crucial features of 
candidate theories of quantum gravity, e.g. via space- 
space uncertainty relations space-time uncertainty 
relations 0, Q-deformed structures ^'^d minimal 

length effects Analyses based on an effective field 
theory approach have also been performed (see 
for a review and a more complete list of references). 

In this paper we work under the hypothesis that super- 
string theory is the correct description of physics on the 
smallest length scales. A key feature which distinguishes 
string theory from point particle-based quantum field 
theory is the new symmetry of T-duality [iJI ■ We will ex- 
plain this symmetry in the context of a theory with closed 
strings only, although the symmetry extends to theories 
with open strings 0) Q ■ Let us assume that each spa- 
tial direction is a torus with radius R. Closed strings have 
various degrees of freedom: momentum modes which cor- 
respond to the center of mass motion of the strings and 
whose energies are quantized in units of 1/i? (setting the 
string length equal to 1 for this discussion), oscillatory 
modes which correspond to wiggles of the strings and 
whose energies are independent of R, and winding modes 
which count the number of times the string winds the 
torus, whose energies are quantized in units of R. 

The string mass spectrum is 



= — + v?R^ + 2[N^ + N^ 



(1) 
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where n is the integer which determines the Kaluza-Klein 
momentum modes, w is the string winding number, and 
N-^ and are integers which denote the number of 
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excited oscillation modes on a closed string in the right- 
moving and left-moving directions around the string. 
Note that in the case of several compact dimensions, n 
and w are vectors of integers, and the squares and w'^ 
appearing above are scalar products of the respective vec- 
tors. The quantum numbers are constrained by the level 
matching condition N-^ — N-^ = nw. In what follows we 
ignore the last term which is only a constant when only 
the leading terms with n,w = 1 are considered. 

T-duality is the transformation which takes the com- 
pactified radius to its inverse, and exchanges the quan- 
tum numbers n and w: 



R 



1 

R 



w . 



(2) 



The mass spectrum is invariant under this transforma- 
tion. The T-duality symmetry has been used to ar- 
gue that a stringy early universe will be nonsingular and 
the special effects of winding modes could also explain 
why there are only three large spatial dimensions (see 
also ^3 for an early paper, and 0, |23| for more re- 
cent developments). This approach to the stringy early 
universe is called "string gas cosmology" . 

To date, effects of T-duality on the spectrum of grav- 
itational waves have not been studied. It is the aim of 
this paper to perform this study in the context of the 
hypothesis that there was a period of cosmological infla- 
tion. Note that at the moment one of the outstanding 
challenges in the field of stri ng gas cosmology is to ob- 
tain a period of inflation (see |2lL 122. l23l | for some recent 
ideas). 

In what follows we propose to study the effects of T- 
duality on the power spectrum of a test scalar field on an 
expanding background (the equation of such a test scalar 
field is equivalent to the equation satisfied by the polar- 
ization modes of gravitational waves) . We will begin with 
a modified propagator which encodes the consequences of 
T-duality |2J| and Fourier transform to obtain a modi- 
fied dispersion relation. We then use the techniques of 
0, 113 to follow the perturbations from the initial time 
(when we assume that they start out in the state that 
minimizes the Hamiltonian) until they cross the Hubble 
radius._ From then on, their evolution is standard (see 
e.g. 26] for a review of the theory of cosmological per- 
turbations and (221 for a pedagogical overview). 



II. PRELIMINARIES: A NEW MODIFIED 
DISPERSION RELATION 

Our starting point is the result of [l^, which shows 
that the T-duality symmetry results in the following 
propagator of the string center of mass (after Fourier 
transforming) : 



1 lsKi{lsy/W 



(3) 



where Kn{z) is the modified Bessel function of order n, 
and Is is the length scale below which trans-Planckian 
physics is important (in our case it is the string length). 
This is the same result found earlier in '28] under the as- 
sumption that quantum gravity effects lead to a modifi- 
cation of the space-time interval on short distance scales, 
the modified formula for the interval being (a; — y)^ + 
instead of {x — y)'^ . 

Making use of the asymptotic forms of the Bessel func- 
tions, the propagator has the following limiting forms 



1 



Greg {k ) 



fc2 + 



(4) 



exp(-ZsVfc^ + m?) , IsV'k'^ -I- m2 > 1 , 



As is evident, for momenta small compared to l^^, the 
propagator approaches the usual one, whereas for large 
wavenumbers it decays exponentially with k. This de- 
cay is a reflection of the fact that T-duality smoothes 
out the usual ultraviolet divergences of the theory. Writ- 
ing the propagator in momentum space as oj"^, where uj 
is the effective frequency, we can read off from Q the 
modified dispersion relation. In the high energy regime 
lsVk^ + rri2 ^ 1 it takes the form 



F{k) = iexp(ZsVfc2 + m2) 



(5) 



Note that in the above (|3J|, I^J and lO, k and w denote 
the physical momenta and frequencies and not the co- 
moving quantities. The comoving wavenumber will be 
denoted by n. 

In cosmology we follow the evolution of fluctua- 
tions which correspond to waves with a fixed comoving 
wavenumber. For these waves, the dispersion relation is 
time dependent. Instead of the simple linear form 



fc2 



it takes the form 



LO = F 



a(ry) 



(6) 



(7) 



where rj denotes conformal time, in terms of which the 
background metric is 



(8) 



dx^ denoting the metric of flat Euclidean three space, 
and 0(77) being the scale factor of the universe. 

To study the effects this modified dispersion relation 
has on the power spectrum of gravitational waves, we 
follow the approach of . We start with a free massless 
scalar field (t>(ri,x.) living in the above background space- 
time. The Fourier modes evolve independently, as do the 
Fourier modes of the cosmological perturbations in linear 
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theory. After introducing the re-scaled field ^ defined via 



(^(?7,x) 



(9) 



(27r)3/2 

the equation of motion for the Fourier modes of reads 



is a more general function of the cosmological background 
(which, in fact, is proportional to a{ri) if the equation of 
state of the background does not change in time - see e.g. 
I2II23 for reviews). 

To take into account the modification of the dispersion 
relation, we need to replace ((Tn|l by 



a 

a 



= 0. 



(10) 



a 
a 



M = 0, 



(11) 



where the effective comoving frequency n^f / is given by 
nlfj = a\v)F'[n/a{rj)] (12) 



Note that gravitational waves in an expanding space obey 
precisely this equation, whereas the scalar metric fluctu- 
ations obey an equation which is very similar. The only 
change is that the time-dependent negative square mass 
— ^ in the above equation is replaced by — where 2(77) 

I 

A. Trans-Planckian Region 

In an inflationary background, any Fourier mode passes through three different intervals of time which are depicted 
in Figure 1. These regions are: 



(RI) Region I: when the physical wavelength of a given 
mode, X{ri) = {2TT / n)a{rj) is smaller than the cutoff 
length where the dispersion relation begins to de- 
viate from its usual linear form. Here, it is taken 
to be the string length Ig'. X < Is- 

(RII) Region II: when the wavelength is longer than Ig 
but smaller than the Hubble radius: Is < X < Ih 

(RIII) Region III: when the wavelength exceeds Ih, the 
mode has frozen out and is being squeezed as a con- 
sequence of the negative square mass term in the 
equation of motion: (i.e crosses the Hubble radius): 
Ih < X. 

A special feature of our dispersion relation is that, in 
Region I (RI) , the non-linearity of the dispersion relation 
leads to a violation of the adiabaticity (WKB) condition. 
This can be quantified in the following manner. We need 
to compare the rate of change of the physical frequency 
i^phys, namely -£^{lnwphys), with the Hubble expansion 
rate Ti = a' /a. If the former is larger than the latter, then 
adiabaticity is violated. Thus, one defines an adiabaticity 
coefficient by: 



a{n,T]) 



n 



cLt] 



(13) 



If a <C 1 then the adiabaticity condition is violated. The 
solution to the modified mode equation has no time to ad- 
just itself such that an initial lowest energy state always 
tracks the instantaneous lowest energy state of the time- 
dependent Hamiltonian of the mode. Thus, even if the 
initial mode wave function minimizes the energy density, 
it is possible that it is in an excited state when it enters 



the second region (RII). In this case, the initial condi- 
tions for the evolution of the mode in Regions II and III 
will be different than what is assumed in the usual linear 
theory of gravitational waves. Since short wavelength 
modes spend more time in the region of non-adiabatic 
evolution than longer wavelength ones, one could expect 
a change in the slope of the power spectrum of cosmo- 
logical fluctuations, as initially argued in |^]. 
For both power law inflation when 



a(?7) 



^0 



l+a 



(14) 



with some constant a > 0, and for exponential inflation 
when 



0(77) 



we see that 



a' , ,1 
n = - = {l + a) — 
a |7y| 



(15) 



(16) 



From Eqs. Q and iQ, and neglecting the mass term for 
now, it follows immediately that 



phys 



= Hklg 



^■phys 

Thus, the adiabaticity coefhcient is 

X_ 

Ts' 



ain,T]) 



(17) 



(18) 



where A is the physical wavelength corresponding to k. 
It follows that the adiabaticity condition is violated in 
Region I and we should hence expect modifications to 
the power spectrum. 
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FIG. 1: Sketch of the relevant regions for our model including the Trans-Planckian region where the WKB approx- 
imation for the mode equation breaks down. The horizontal axis is physical distance, the vertical axis conformal 
time. The red (dotted) line indicates the Planck or string length, the distance where the usual dispersion relation 
becomes invalid. The green (solid) line which is vertical during the period of inflation which lasts from rji to rjn, 
is the boundary between Regions II and III, the regions where the mode evolution is oscillatory or frozen, respec- 
tively. The two thin lines for which the distance grows exponentially during the period of inflation indicate the 
physical wavelengths of two fixed comoving modes fci and /c2. As is evident, they start out in Region I, the Trans- 
Planckian region, at early times, and then cross Region II to end up in Region III during the later stages of inflation. 



B. Minimum Energy Density Initial Conditions sity of the field. With this prescription, the initial value 

of the mode function becomes 

InRcf.1, following Ref. |29^ . it was argued that since 
the WKB approximation is not valid in the Trans- 
Planckian regime, we cannot choose as our initial con- 
ditions the Bunch-Davies vacuum. Instead we have to _ ^ -j _ 1 _ ^ ^^g^ 
choose the state which initially minimizes the energy den- ^j2Lo{rji) -y/SnT// 
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and from the minimal energy density requirement we ob- 
tain: 



fJ-'iv = Vi) = ±* 



(20) 



Here, rji is the initial time at which we choose to minimize 
the energy density. 

Other relevant mode-dependent times are: 

(1) ryi, the time when the mode labeled by a fixed co- 
moving wavenumber n crosses from Region I to Re- 
gion II, is given by 



l+a _ k 27r 



Is n 

thus defining a constant e ^ 1 



27r 
ne 



(2) ri2, the time when the mode crosses the Hubble 
radius, i.e. exits from Region II into Region III, is 
given by 



27r 



\m\ 



(1-1- a) 



(22) 



In Region I, the argument in the exponential in (|23|l is 
large, i.e. 2ne|77| > 1, and hence we can approximate the 
equation of motion in the following way: we can drop the 
second term - the squeezing term - in the square paren- 
theses of (I23|) since it is suppressed compared to the first 
term. In the first term, we replace rj by 771 in the denomi- 
nator since the time dependence in the exponential over- 
whelms the time dependence of the other factors. Thus, 
the approximate form of the equation of motion Eq. 124() 
in Region I becomes: 



o2ne|i)| 2 



47r2 



-^i{r,) = 0. 



(21) After the change of variables 



exp(rie|77|) 
2^e 



the equation (|24|l takes the form 



x'^li" + a;/i' -I- a;^/i = , 



(24) 



(25) 



(26) 



III. SPECTRUM IN THE CASE OF MINIMAL 
ENERGY INITIAL CONDITIONS 



which is nothing but Bessel's differential equation of ze- 
roth order. The solution can be written as 



In this section, we present approximate analytical solu- 
tions for the mode functions obtained using the modified 
dispersion relation derived above, and compute the re- 
sulting power spectrum. For each mode, we solve the 
mode equation in each of the three regions (see Figure 
1) and apply the standard matching conditions at the 
boundaries between the regions. 



A. Region I (RI) 

Inserting the dispersion relation H12|l into the mode 
equation and making use of jSJ and of the form of the 
scale factor for an exponentially expanding background 
p5|l - for simplicity we will focus on this case, i.e. on 
a = - we obtain 



exp(2?ie|77| 



a 
a 



(23) 



(27) 



where Jo and Iq denote the Bessel functions and A and 
B are constants. 

Since we are interested in the high energy behavior of 
this solution (recall that we are in the Trans-Planckian 
regime) , we need the asymptotic representations for large 
arguments of the zeroth-order Bessel functions. 



cos(a; - ^)J , V X > 1 (28) 



Joix) ~ W — 



/ 2 r , 

Yo{x) ~ W — sm a; - - 

V TTZ L 4 



, Va;> 1. 



(29) 



Recalling the expression for x, the asymptotic solution 
/i = /i/ in (RI) reads 



Ai cos 



27re 



7r/4 



A-} sin 



27re 



7r/4 



(30) 



Note that in terms of the variable x, the two distinguished mode-dependent times 771 and 772 take the form 



= = (31) 

27re 27re 



gn«l'72| 
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In terms of x, the solution reads 
"2 



(Ai co8{x - 7r/4) + A2 sm{x - 7r/4)) 



(33) 



conditions. To implement these, we compute the deriva- 
tive of fi and use 

d dx d d 



The coefficients Ai and A2 are determined by the ini- 
tial conditions. We choose minimal energy density initial From Eqs. H28|l and (|29|l 



and 



rj(x) ^ ne 



sin(x — 7r/4) H cos(a; — 7r/4) 

2x 



cos{x — 7r/4) — — sin(x — 7r/4) 



and the asymptotic form of /i^ becomes 



n'j{x) ~ neA / — [—^1 sin(a: — 7r/4) + A2 cos{x — 7r/4)] 



(35) 



(36) 



(37) 



To solve for Ai and A2 we need to apply the minimum energy density conditions Eqs. H19f) and (|20|l . but cast in 
terms of Xi instead of |r/i|. 



Ai cos{xi — n/A) + A2 sin(xi — 7r/4) 



-Ai sin(a;i — 7r/4) + A2 cos(xi — 7r/4) = ±- 



nxi 



An. 



ne V 4x, 



(38) 
(39) 



We can now use the fact that e <C 1 and the condition that ne\ri\ ^ 1 to guide us through our approximations. Solving 
for Ai and A2 we get 



Ai 

A2 



cos(a;i — 7r/4) =F — . / — sm(xi — 7r/4) 



An, 



ne V Axi 



sm(Xi — n/A) ± — . coa{xi — tt/A) 



Aneffi ne V Ax 
Comparing the two terms in each expression, we find that 

Ai 

A2 ~ J , sin(xt - 7r/4) 



An 



eff. 



■ cos{xi — n/A) 



An 



Eq. (123 thus gives us the following solution in Region I 
fj.i{x) 



1 Xi 

— \ — [cos(xi — n/A) cos{x — n/A) + sm{xi — n/A) sm{x — n/A)] 

V^neff^ V X 



(40) 
(41) 

(42) 
(43) 

(44) 



or, more compactly, 

Hi{x) = 



— cos(xi — x) 



B. Region II (RII) 

(45) In Region II, the adiabaticity condition is not violated, 

and thus the mode functions /i = /i// are oscillatory, i.e. 
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by 



^lu{^l) = Bie"" + Bse"'"" (46) 
or equivalently in terms of x 



fi'jj{x) 



Bi{2Trey/'x'/' - B2{2ney 



i/e 



X 



.(48) 



The coefScients Bi and B2 are determined by the 
matching conditions at the time \r]i \ = ^ 



f^i{\m\) = f^n{\Vi\) 
/^/(hi|)-/^/7(hi|) 



(49) 



or their equivalent in terms of the variable x. Recall that, 
in terms of the x variable, fii{xi) and fi'j{xi) are given 



^^I{xl) = ^J^^^^x-^^^ coa{xi - x^) (50) 

fl'j{x) ~ 



ne In \ 2TTeXj 



Air 



sm{xi — Xi) . 



(51) 



where we have again used the fact that e <C 1 and ne\ri\ ^ 
1. To abbreviate the notation, wc introduce the variable 



Ai = y/hi\2nexi 



(52) 



Then, the matching conditions (|49|l (making use of the 
explicit value of xi from Ea. (|^ ) read 



V2n 



A^ cob{x, - xi) = Bi{2TTeyl^xl^ + B2{2iiey^x^ 



^/2J^ 27r 

The coefficients Bi and B2 are thus given by 



sin(a;, - xi) = Si(27re)*/^a;*/^ - B2{2TTt)-'''x-'/' 



Bi = 



B2 = 



2V2n 



2J2n 



{{2nerl^xT'' 
({2Tiey''xl' 



e cos(xi — xi) — i— sin(a::i ~ Xi) 
2tt 

e ^ cos{xi — xi) + i— sin(xi — xi) 
2tt 



(53) 
(54) 

(55) 
(56) 



This time, we must keep all the terms in the coefficients since they are all of the same order. Thus, the solution of 
the mode function in (RII), Ea. H47() . reads 



fiii{x) 



2V2n 
or, more compactly. 



X y'" f ^ 

— I [ G ^ coslxi — xi) — i — sm(xi — Xi 
xi J \ 2tt 



e ^ cos(xi — a;i) + i— sin(a;i — Xi] 
2tt 



A 



e cos(xi — xi) cos 



In \x/xi I 



e'' . /ln|a;/a;i| 

— sm(xi - xi)sm 

Ztt \ e 



(57) 



r 



C. The Power Spectrum 



Now we are in a position to compute the power spec- 
trum. Assuming 0(1) coupling between the growing 
mode in Region III and /1/7, the non-decaying mode mode 
in Region III is 



Hiiiir)) ~ Ca{r)) 



(58) 



where the constant C is fixed by matching /i/// with /i// 
at the Hubble-crossing time I772I (see Ea. (|^ ') 

lo 



f^iii{\m\) - Ca{\rj2\) = C- 



(59) 



In terms of the variable x the matching conditions are 



and hence 



tJ-iiiix2) = Cn^ = M//(a;2), 
Ztt 



C = —rtJ'ii[X2) ■ 

nlo 



(60) 
(61) 



8 



The power spectrum for our scalar field is defined Making use of (|61|l and (|57|) we get 



as 



P^in) = n3|(/)(fc)P = n3a-2|;,(n)|2. (62) 

Tfie spectral index tit for gravitational waves is defined 
by 



where At is the amplitude of the spectrum. 
Combining Ht)2|) and (|58|l we obtain 



(63) 



P^{n) = Mi^|^,,(|^2|)|' (65) 

= -^\AA^J{\m\Mi\M2\)? (66) 



P^{n) = 



(64) 



where i/dry^l, \t]2\) will yield oscillations in the spectrum of a particular form, and is given by 



J{\-ni\, \m\) = e ''cos 



27re 



cos 



ne\rii\ 



271 



■ sm 



g«e|')i| _ gn«l'7l| 



27re 



sm 



ne\'q2\ 
ne\'qi\ 



(67) 



Plugging in the values for \rii\ and |?72| and making more transparent the ne|ryi|-dependence of j/(n, e, |?yi|), we obtain 



^"(71, e, Iry^l) = cos 



' sm 



27re 



27re 



e cos 



27rey 27r 

„27r 



2tt€ 

2-K - 



— 1 sm cos 

\2TifJ 2ti \2TTel 



(68) 



The amplitude and spectral index of the power spec- 
trum can be read off by inserting the expression (|52() for 
Ai and making use of the definition (|25|l of the variable 
X, yielding 

Axn"- = ^ ne|77,| \ J{n, e, \r^,\)\^ . (69) 



of matter on scales of the order of the atomic separa- 
tion resulting in a departure from the linear dispersion 
relation. 



D. Discussion 

From our main result (gni) we draw the following con- 
clusions. Given our modified dispersion relation for fluc- 
tuations, the spectrum of gravitational waves in an ex- 
ponentially expanding background is characterized by a 
non- vanishing spectral tilt of nx = 1. In addition, there 
are oscillations in the spectrum with a characteristic de- 
pendence on n. This result was obtained working to lead- 
ing order in e. 

The fact that the spectrum is not scale-invariant but 
blue-shifted is due to the fact that the shorter wave- 
lengths are subject to the modified dispersion relation 
for a longer time, leading to a higher excitation level of 
the modes when their wavelength is equal to the cutoff 
scale after which the dispersion relation becomes linear. 
This is similar to what happens in models with a dis- 
crete space-time, or in analog models from condensed 
matter physics where the modes feel the granular nature 



We would like to point to an important observation 
regarding Ea. (|59|l . We made it explicit to show how the 
combination ne|7/i| appears both in the oscillatory func- 
tion J' and in Eq. (|69() . We see that to get a scale invari- 
ant spectrum without oscillations, it suffices to modify 
the prescription for the initial conditions and postulate 
that the modes are created at a mode-dependent initial 
time \r]i \ for which n\r]i \ = constant (this is similar to what 
is postulated to occur in the analyses of 0? ]). We will 
discuss later on a better derivation of what the initial 
conformal time might be by linking the fact that hav- 
ing T-duality at the basis of our dispersion relation, i.e. 
as our Trans-Planckian physics, is similar as modifying 
our space-time metric by adding to it a minimal length 
which is of the order of the string length Is- This could 
be linked to non- commutative space-time or to a stringy 
space-time uncertainty relation which was studied earlier 
in Ref . P| . We will return to that point at the end of our 
article since it is of a more speculative nature. 
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IV. SPECTRUM IN THE CASE OF AN 
INSTANTANEOUS MINKOWSKI VACUUM 

To test the dependence of our previous result on the 
initial conditions we adopt another choice of initial con- 
ditions, the instantaneous Minkowski vacuum as an al- 
ternative to Eas. lfTt^ and f^. That is, at \ri\ = \rii\ the 
mode function now satisfy: 



1 



= ± 



2n 

n 
2 



(70) 



(71) 



A. The Trans-Planckian Region (RI) 

Following the same steps as in the case of the minimum 
energy density initial state, we set the initial conditions 
at a fixed time rji in the Trans-Planckian region (RI) . We 
start with the general solution in Region I from Eas. (|33(l 
and ()37|l and match them to the initial values H70|l in 
order to determine the new coefficients Ai and A2. The 
result is 



^T- - y cos{xi - 7r/4) 

, / IT X 2 . / I A \ 

A2 ^ W -— sm(a;i - 7r/4) , 
V 4n 

and the mode fimction fJ.i{x) is thus given by 
1 



fii{x) 




cos{xi — x) 



(72) 
(73) 

(74) 



B. Region II 

In Region II where the WKB approximation is valid, 
we have plane wave solutions for the mode functions as 
in Eqs. H47|) and H48|l. The mode functions are matched 
at the time |77i| when the physical wavelength of a mode 
equals the minimal length Is, that is 



fj^iixi) 
n'lixi) 



fJ-iiixi) 
f^'iiixi) 



(75) 
(76) 



These yield two equations for the coefficients Bi and B2 
which can be solved to get 



^ ~i jxiXi sin(x 
^ ~ Y\l^ (27rexi)*/ 



Xi) 



Bo 



2n 



sm{xi — Xi){2'KexiY^'^ 



(77) 
(78) 



The mode solution fijj (x) in Region II is thus given by 

In 



fj,ii{x) = 



2n 



- sin(a;i — xi) sin 



(79) 



C. The Power Spectrum Revisited 

Again assuming that the growing mode in Region III 
picks up 0(1) of the amplitude of nu at the matching 
time X2, the solution /i/// in Region III is given by 



fiiii = Ca{r]) , 



(80) 



with the constant C determined by C = ]^/i//(x2) or, 
expanded out 



2tt XiXi 
C = - — \l— — sm(a;i — a;i) sm 



ln|x2/a;i| 



Ion V 2n 

The power spectrum P,j,{n) thus becomes 
P^{n) = Asn"^ = n^lCp 

27r2 2. N 2 

'0 



In terms of conformal time, 



-2lf- 



27re 



ln|a:2/xi| 



27r 



(81) 



(82) 

■ :83) 



(84) 



D. Discussion 



The striking difference between our previous result ob- 
tained using the minimum energy density calculations 
and that found here using the instantaneous Minkowski 
vacuum initial conditions is the dependence of the power 
spectrum on the comoving wave number n. We see from 
Eq. (|84|l that the dependence on n is exponential while 
from Eq. (|69|) the dependence on n is linear. Thus, for 
the instantaneous Minkowski vacuum we get an exponen- 
tially blue spectrum (note that the quantity ne|ryi| which 
appears in the exponential is much greater than one). 

If, instead of imposing initial conditions for all modes 
at a fixed time, we use the alternative discussed section 
nil Dl namely assume a minimal length in our theory and 
argue that initial conditions should be set on a surface of 
constant n\rii\ which would correspond to the surface on 
which the physical mode are created, and setting 



n\r]i 



ln|27r2e(5| 



(85) 



(for some fixed number S) we would obtain the following 
spectrum: 



o2-7r^2 



■ sm 



tt5 - 



27re 



2-K 



(86) 



which is scale invariant spectrum without any oscillations 
which depend on the comoving wavenumber n. Thus, for 
this initial time prescription we arrive at the same conclu- 
sion as in the case of minimal energy initial conditions. 
This is not a real surprise since in this case all modes 
spend the same amount of time in the trans-Planckian 
region. 
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V. CONCLUSIONS 

In this paper we have studied the spectrum of gravi- 
tational waves which results when the usual linear dis- 
persion relation for the wave equation is replaced by the 
non-linear dispersion relation which results when assum- 
ing that the wave propagator is consistent with the T- 
duality symmetry of string theory. This modified disper- 
sion relation differs from the linear one on length scales 
smaller than the string scale. 

We have shown that the modified dispersion relation 
leads to a non-adiabatic evolution of the mode functions 
in the trans-Planckian region. If, as is usual in an in- 
flationary cosmology, the initial conditions for fluctua- 
tions are set for all modes at the same initial time, this 
modification, in turn, leads to a change in the slope of 
the spectrum of gravitational waves. Instead of a scale- 
invariant spectrum, a spectrum with blue tilt given by a 
spectral index nx = 1 results. In addition, the spectrum 
has characteristic oscillations. 

Our results are based on studying the evolution of 
a test scalar field on our expanding inflationary back- 
ground. According to the theory of cosmological pertur- 
bations, each gravitational wave polarization mode obeys 
the same equation as such a test scalar field. The scalar 
metric fluctuations, in turn, obey an equation with a 
slightly different squeezing factor [23, |23| ■ The difference 
in the equations of motion, however, is in general only 
important on scales larger than the Hubble radius, and 
even then only if the equation of state of the background 



is changing in time. Thus, it is reasonable to expect that 
a similar analysis to the one given in this paper applies 
to scalar metric fluctuations. The result for scalar met- 
ric fluctuations would then be a spectrum with blue tilt 
ns — 2, inconsistent with recent data js^l- If the only 
change to the equation of motion for cosmological fluctu- 
ations in the context of string theory on trans-Planckian 
scales were the change in the dispersion relation discussed 
in this paper, our results would imply a deep problem in 
realizing successful inflation in the context of string the- 
ory. However, as discussed in |||, there are other aspects 
of string theory, e.g. the space-time uncertainty relation, 
which must be considered, and these considerations could 
well restore the prediction of a scale-invariant spectrum 
in the case when the background space in exponentially 
expanding. 

In conclusion, we hope to have convinced the reader 
that, assuming that there was a period of cosmological 
inflation, the basic principles of string theory are clearly 
testable in cosmological observations. Final predictions, 
however, will have to await a better understanding of the 
equations of string theory on trans-Planckian scales. 
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